Sodium diffusion through amorphous silica surfaces: A molecular dynamics study 
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(Dated:) 

We have studied the diffusion inside the silica network of sodium atoms initially located outside 
the surfaces of an amorphous silica film. We have focused our attention on structural and dynamical 
quantities, and we have found that the local environment of the sodium atoms is close to the local 
environment of the sodium atoms inside bulk sodo-silicate glasses obtained by quench. This is in 
agreement with recent experimental results. 
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I. INTRODUCTION 

Many characteristics of materials such as mechani- 
cal resistance, adsorption, corrosion or surface diffu- 
sion depend on the physico-chemical properties of the 
surface. Thus, the interactions between the surfaces 
with their physico-chemical environment are very impor- 
tant, and in particular for amorphous materials which 
are of great interest for a wide range of industrial and 
technological applications (optical fibers coating, cataly- 
sis, chromatography or microelectronics). Therefore a 
great number of studies have for example focused on 
the interactions between the amorphous silica surfaces 
with water, experimentallyi and by molecular dynamics 
simulationo 2 ! 3 ! 4 ! 5 ! 6 ! 7 . 

On the other hand, the sodium silicate glasses entail 
great interest due to their presence in most of the com- 
mercial glasses and geological magmas. They are also 
often used as simple models for a great number of sil- 
icate glasses with more complicated composition. The 
influence of sodium atoms on the amorphous silica net- 
work is the subject of numerous experimental studies: 
Raman spectroscopy^, IR&ii, XPS±2J£ and NMRiiii 
from which we have informations about neighboring dis- 
tances, bond angle distributions or concentration of so- 
called Q™ tetrahedra. In order to improve the insight 
about the sodium silicate glass structure and to obtain 
a good understanding of the role of the modifying Na + 
cations, Greaves et al. have used new promising investi- 
gation techniques like EXAFS and MAS NMRi£J&. De- 
spite all these efforts, the structure of sodo-silicate glasses 
is still a subject of debate. Another means to give infor- 
mations about this structure is provided by simulations, 
by either ab initial! or classicali2ii2i22i2ii22i2ii molecular 
dynamics (MD). In the present work, we are using clas- 
sical MD simulations, but a contrario to previous sim- 
ulations, the sodium atoms are not located before hand 
inside the amorphous silica sample. 

Recent experimental studies of the diffusion of Na 
atoms initially placed at the surface of amorphous silica, 
using EXAFS spectroscopy^, showed that the Na atoms 
diffuse inside the vitreous silica and once inside the amor- 



phous silica network, the local environment of the Na 
atoms is characterized by a Na - O distance G?Na-o = 2.3 
A and by a Na - Si distance G?Na-Si = 3.8 A. These values 
are close to the distances characterizing the local environ- 
ment of Na atoms in sodium silicate glasses obtained by 
quench. 

In this study we have used classical MD simulations in 
order to reproduce the diffusion of sodium atoms inside a 
silica matrix and to check that the local environment of 
the sodium atoms is close to what is found for quenched 
sodo-silicate glasses. The sodium atoms have been in- 
serted at the surface of thin amorphous silica films under 
the form of Na20 groups in order to respect the charge 
neutrality. 



II. COMPUTATIONAL METHOD 

To simulate the interactions between the different 
atoms, we use a generalized version^ of the so-called 
BKS potential— where the functional form of the poten- 
tial remains unchanged: 
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The potential parameters Aij , Bij , dj , qt and qj involv- 
ing the silicon and oxygen atoms (describing the interac- 
tions inside the amorphous silica network) are extracted 
from van Beest et al. 2e and remain unchanged (in partic- 
ular the partial charges qsi = 2.4e and qo = — 1.2e are 
not modified). The new parameters, devoted to describe 
the interactions between the sodium atoms and the silica 
network are given by Kramer et alm^ and are adjusted 
on ab initio calculations of zeolithes except the partial 
charge of the sodium atoms whose value qNa = 0.6e is 
chosen in order to respect the system electroneutrality. 
However, this sodium partial charge does not reproduce 
the short-range forces and to this purpose, Horbach et 
alm& have proposed to vary the charge qNa as follows: 
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where r,j is the distance between the particles i and j. 
The parameters C and r c are adjusted to obtain the ex- 
perimental structure factor of Na 2 Si 2 0s (NS2) and their 
values are included in Ref 22 . It is important to note that 
using this method to model the sodium charge, the sys- 
tem electroneutrality is respected for large distances (in 
fact for distances r ^ r c ). Next we assume that the mod- 
ified BKS potential describes reasonably well the system 
studied here, for which the sodium atoms are initially lo- 
cated outside the amorphous silica sample. In addition 
other simulations have shown that this interatomic po- 
tential is convenient for various compositions, in partic- 
ular for NS2, NS3 (Na 2 Si 3 7 )^ and NS4 (Na 2 Si 4 9 )S 
and we assume it is adapted to model any concentration 
of modifying Na + cations inside sodo-silicate glasses. 

Our aim here is to obtain a sodo-silicate glass by depo- 
sition of sodium atoms at the amorphous silica surface, as 
it was done experimentally 24 . Using the modus operandi 
described in a previous study2& we have generated Amor- 
phous Silica Films (ASF), each containing two free sur- 
faces perpendicular to the z-direction. These samples 
have been made by breaking the periodic boundary con- 
ditions along the z-direction, normal to the surface, thus 
creating two free surfaces located at L/2 and — L/2 with 
L = 35.8 A. In order to evaluate the Coulomb interac- 
tions, we used a two-dimensional technique based on a 
modified Ewald summation to take into account the loss 
of periodicity in the z-direction. For further technical de- 
tails see Re£2&. Then, instead of initially positioning the 
sodium atoms inside the silica matrix, like it was done 
beforeiSa^lLSiiSiSi, we have deposited 50 Na 2 groups in- 
side two layers located at a distance of 4 A of each free 
surface as depicted in FIG.Q] 
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FIG. 1: Initial configuration: two layers containing 50 Na20 
groups each placed at 4 A of the two free surfaces of an amor- 
phous silica film. The way the Na2 groups are arranged 
within the layers is shown in the blow up. 

Within the layers, the Na 2 groups are assumed to be 



linear, with dwn-o = 2 A2&, and arranged on a pseudope- 
riodic lattice represented in the zoom of FIG. Hence 
the system is made of 100 Na 2 groups for 1000 Si0 2 
molecules, corresponding to a sodo-silicate glass of com- 
position NS10 (Na 2 Siio0 2 i) and contains 3300 particles. 
Since our goal is to study the diffusion of the sodium 
atoms placed at the amorphous silica surfaces, we fixed 
the initial temperature of the whole system at 2000 K. 
Indeed, the simulations of Smith et alm& and Oviedo et 
alw^ of sodo-silicate glasses have shown that there is no 
appreciable sodium diffusion for temperatures below ss 
1500 K. On the other hand, it is worth noticing that 
Sunyer et alM> have found a simulated glass transition 
temperature T g ~ 2400 K for a NS4 glass. Therefore 
we thermalized the sodium layers at 2000 K and placed 
them at the ASF surfaces, also thermalized at 2000 K. 

We have used repulsive 'walls' at z = —30 A and 
z = +30 A in order to avoid that some Na 2 groups 
evaporate along the z direction, where the periodic 
boundary conditions are no more fulfilled. The en- 
ergy of the repulsive 'walls' has an exponential shape, 
E = Eq cxp [— [z w — z) /a) where z w is the wall position, 
(7 = 0.1 A the distance for which the repulsion energy is 
diminished by a factor e and Eq = 10 eV the repulsion 
energy at the plane z = z w . The value of 30 A for z w was 
chosen in order to place the repulsive walls at a reason- 
able distance from the Na 2 layers and not too far from 
the ASF surfaces. 

We have then performed classical MD simulations, 
with a timestep At = 0.7 fs, using ten statistically in- 
dependent samples. Since the interactions between the 
surfaces and the Na 2 layers are relatively weak, some 
Na 2 groups may evaporate just before being reflected 
toward the thin films by the repulsive walls. During this 
time frame, the system temperature increases up to a 
temperature of approximately 2800 K. As described by 
Athanasopoulos et al£^ this temperature rise is likely 
due to the approach of the adatoms to the surface of 
the substrate, dropping in the potential well of the sub- 
strate atoms, thus increasing their kinetic energy. This 
kinetic energy is then transmitted to the substrate with 
the adsorption. Contrarily to Athanasopoulos et al. 31 
and Webb et al£2> we have not dissipated this energy 
excess with a thermal sink region, in fact we have not 
controlled the temperature at all. Using this device we 
are able to perform MD simulations of the diffusion of the 
sodium atoms deposited via Na 2 groups at the ASF sur- 
faces. In the following section, we will present the struc- 
tural and dynamical characteristics of the sodo-silicate 
film (SSF) obtained in this way. 



III. RESULTS 

The observation of the sodium diffusion in the silica 
network is an important goal of this molecular dynam- 
ics simulation. This can be carried out by analyzing the 
behavior of the density profiles along the normal direc- 
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tion to the surface (the z-direction). The density pro- 
files represent the mass densities within slices, of thick- 
ness Az = 0.224 A, parallel to the surfaces 28 . The time 
evolution of the sodium density profile is represented in 
FIG. Ufa) and FIG. H» and the time evolution of the 
total density profile in FIG. Etc) and FIG. Eld) . 
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FIG. 2: Time evolution of the Sodium (a) and (b) and total 
(c) and (d) density profile. 

The sodium atoms enter inside the silica network in the 
time range 14 - 70 ps (FIG. Eta)) and during this time in- 
terval the sodium profile exhibits a diffusion front. After 
42 ps, the sodium atoms are observed in the entire sys- 
tem illustrating that they have diffused within the whole 
ASF, as observed experimentally^. This result contrasts 



with that of MD simulations of the diffusion of platinum 
(electrically neutral)3i and potassium2£i2LSi but is sim- 
ilar to that of the diffusion of lithiurnSiiSLSi at the sur- 
face of amorphous silica. Moreover, the sodium profiles 
(FIG.Et a ) and (b)) show outward rearrangement at the 
surface as already pointed out by Zirl et alM for glassy 
sodium aluminosilicate surfaces. This agrees with ion 
scattering spectroscopy 3 ^ and simulations of sodium sil- 
icate glasses 3 ^ in which high concentrations of sodium 
ions are found at the surface. For larger times, t ^ 98 ps 
represented in FIG. Etb), the sodium density seems to 
stabilize around a mean value of 0.3 g.cm -3 which cor- 
responds to ~ 2 Na atoms per slice (Az = 0.224 A ). 

The surface of the system is identified as being the 
large linear region in which the total density profile de- 
creases. For the short times (FIG. Et c )) the surface is 
located approximately in the range 15 - 20 A and for 
larger times (t ^ 98 ps, FIG. Eld)), the surface lies in the 
range 7 - 25 A. Hence, the introduction of the sodium 
atoms in the ASF is likely to increase the surface thick- 
ness of the system. On the other hand, as observed for 
the adsorption of platinum atoms on the surface of amor- 
phous silica 3 ^, the surface position does not seem to vary 
with time. We have also calculated the silicon and oxy- 
gen density profiles, but since they behave like the total 
density profile, they are not represented here. The non 
bridging oxygen (NBO) density profile is not represented 
in FIG.[3as well since it is close to the Na density profile. 

In the time range 98 - 210 ps, the total and sodium 
density profiles do not evolve with time. In particular, in 
the region z < 5 A (FIG. Eld)), the total density value 
remains fluctuating around a mean value of 2.6 g.cm -3 
and the atom composition is approximately of 370 Si, 760 
O and 40 Na which is usually written Na 2 (SiC>2)i8.5 or 
'NS18.5' (it should be noticed that the above mentioned 
density is significantly larger than the one expected for a 
"real" NS18.5 glass (w 2.3 g.cm -3 )). Therefore it seems 
reasonable to consider that the system is in a quasi per- 
manent regime after 210 ps and the forthcoming quan- 
tities, structural and dynamical, are calculated for the 
following 70 ps. It is worth remembering that all the 
quantities are determined for a system containing 3300 
particles and averaged over 10 statistically independent 
samples. 

As usual when studying the structural and dynami- 
cal characteristics of free surfaces, the system is divided 
into several subsystems: here six slices of equal thickness 
10 A. But in order to increase the statistics, the con- 
tributions to the physical quantities of the negative and 
positive slices are averaged. Hence, the system is actually 
subdivided into three parts, named respectively from the 
center to the surface, interior, intermediate and external 
region. 

In order to improve the characterization of the local 
environment of the atoms, we have calculated the ra- 
dial pair distribution functions for all the pairs <E 
[Si, Na, O] 2 within the three subsystems defined previ- 
ously. The Na - Na, Na - O and Si - Na pair distribu- 
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tion functions are represented in FIG.|3Ja),[3{b) andEJc) 

respectively and represent the local environment of the 
sodium atoms for t > 210 ps. 

At the surfaces of the system the distances are 
c?Na-o — 2.2 A, dsi-Na — 3.5 A and, due to a lack of 
statistics, the Na - Na distance is included in the interval 
3-3 < dNa-Na ^5 3.9 A. While slightly smaller, these 
distances are close to the experimental values found by 
Mazzara et al. (c^a-D = 2.3 A and dsi-Na = 3.8 A)24. 
Moreover, the values found in the present work agree 
with the distances, calculated by MD, corresponding 
to the sodium environment in sodo-silicate glasses, 
of several sodium compositions (NS3i&*2£, NS3^ and 
NS4ii), obtained by quench. Therefore, as observed 
experimentally by Mazzara et alm^, once within the 
amorphous silica network, the sodium atoms have the 
same local environment as in the sodo-silicate glass 
obtained by quench. This fact is confirmed by the 
distributions of the NaONa and SiONa bond angles 
(not shown) which are close to those determined in 
quenched sodo-silicate glasses^i* 3 ^. In particular, the 
most probable angles are 90° for NaONa and 105° for 
SiONa, in agreement with the values found by Oviedo et 
alm^ and Sunyer et alm&. 
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FIG. 3: (a) Na - Na, (b) Na - O and (c) Si - Na radial pair 
distribution functions for the (• • • ) interior, ( — ) intermediate 
and ( — ) external region. 

The intermolecular distances corresponding to the 
amorphous silica network structure are cZsi— si — 3.1 A, 
dsi-o — 1-6 A and do~o — 2.6 A. On the other hand, 
the most probable values exhibited by the distributions 
of the OSiO and SiOSi bond angles (not shown) are re- 



spectively ~ 109° and ~ 145°. These distances and bond 
angle distributions are very similar to those found ex- 
perimentally or by MD simulations in bulk amorphous 
silica. The shoulder exhibited at 2.5 A by the distribu- 
tion gsi_si a t amorphous silica surfaces 28 and interpreted 
as the signature of the twofold rings is not present in the 
SSF. The absence of the shoulders_jt_8° and_100° de- 
picted at ASF surfacea 2 i 28 i 38 by the OSiO and SiOSi bond 
angle distributions confirms that the small sized rings 
have disappeared as suggested by the radial pair distri- 
butions. This result is expected for R2O (with R belong- 
ing to the column I) adsorption on glassy silicate surfaces 
and observed via MD simulations for R = H, Li, Na and 
K 2 ! 3 ! 6 ! 35 ! 36 ! 37 ! 39 . Note that this also occurs for platinum 
adsorption on sodium aluminosilicate surfaces^i^S and 
for H2 adsorption on amorphous silica surfaces^. The 
small rings like two and threefold rings are known to be 
some of the most reactive sites on the surfaces of silicate 
glasses since they include strained siloxane bonds that 
react with water or other adsorbate o 40 ' 41 . 

Since the sodium introduction weakens the amorphous 
silica network, one important question is to measure the 
proportion of non bridging oxygens (NBOs). To this pur- 
pose, oxygen coordination with silicon was calculated. As 
expected, there is an important proportion of defective 
oxygens due to the presence of modifying cations Na + 
(8.9 % for the SSF to be compared with 1 % for the 
ASF 28 ), illustrated by the similarities between the NBO 
and Na density profiles mentioned previously. Moreover, 
the NBO concentration is coherent with the concentra- 
tion of 11 % of NBO calculated by MD simulations in 
a NS9 system^ 7 .. The latter concentration is higher than 
the one found in this study because of the greater propor- 
tion of Na atoms in NS9 compared to 100 Na 2 for 1000 
Si0 2 in the present SSF. At the surface (z ^ 20 A), the 
BOs disappear and correlatively, the defective oxygens 
become preponderant, as observed for amorphous silica 
(36.2 % of NBOs at the SSF surfaces, 15 % at the ASF 
surfaces^ and 10 % at the nanoporous silica surfaces 4 ^). 
When analyzing the silicon coordination with oxygen, we 
can state that the silicon atoms remain coordinated in a 
tetrahedral way, revealing that the sodium introduction 
does not modify the silicon environment. The modifica- 
tions created by the Na + cations are not able to break the 
Si04 tetrahedra which are very stable in the amorphous 
silica network. This agrees with the usual models for the 
sodo-silicate glasses, ie the CRN model of Zachariasen 43 
and the MRN model of Greaves^. 

One possible way to analyze more precisely the silicon 
environment consists in calculating the Q n tetrahedra 
proportion. A Q™ structure is a Si04 tetrahedra which 
contains n BOs. The Q™ proportion is often determined 
by NMR experiments 44,45 , XPS 13 ^^ or by molecular dy- 
namics simulations^, in order to describe the local en- 
vironment around the silicon atoms. In the SSF, the Q 2 
and Q 3 concentrations (6.8 % and 25.4 % respectively) 
are weak compared to those determined for NS2, NS3 
and NS4 glasses. This result is coherent since the NS2, 



NS3 and NS4 glasses contain more Na + cations than the 
SSF studied in this work. At the surface, the Q 3 pro- 
portion is 45.9 % which is comparable to the proportions 
obtained by MD (using the BKS potential) in NS2, NS3 
and NS4 glasses and to the experimental proportions in 
NS 3iii4 and nsi 4 ^ 4 ? glasses. Moreover, it is worth not- 
ing that some Q 1 appear at the surface. In fact, these 
structural entities do not allow to create a network but 
it is conceivable to find those defects forming 'dead ends' 
at the surface. 

A direct method to confirm the previous assumption 
concerning the disappearance of small rings consists in 
analyzing directly the ring size distribution. A ring is 
a particularly interesting structure because it can be 
detected using infrared and Raman spectroscopy. In 
particular the highly strained twofold rings result in 
infrared-active stretching mode a 41 ' 47 at 888 and 908 
cm . In order to determine the probability P n for 
a given Si atom, whose coordinate along the normal 
direction to the surface is in one of the three regions, to 
be a member of a n-fold ring we have used the algorithm 
described ini 8 .. A ring is defined as the shortest path 
between two oxygen atoms, first neighbors of a given 
silicon atom and made by Si - O bonds. The ring size is 
given by the number of silicon atoms contained in the 
ring. 

The probability P n is reported in FIG. Ufa) for n — 
2, . . . , 9 and for the three different regions. For compar- 
ison, we have also reported P n determined at the ASF 
surface and interior^. In order to improve the medium 
range order characterization, we have investigated the 
orientation of the rings computing < cos 2 > for a given 
ring size, within the three regions of the system, where 
8 is the angle between the normal of the surface and 
the normal of the ring 2 -. The results are reported in 
FIG. Efb) for the three regions and for n = 2,..., 9 
together with the results obtained for the ASF surface 
(dashed line) and interior (dotted line). 

For the three regions, the distributions of the sodo- 
silicate film (FIG. Eta)) are closer to the ASF interior 
than to the ASF surface distributions. Particularly, the 
probability of a silicon atom to belong to a small sized 
ring (2, 3 or 4-fold ring) is weak. This confirms the pre- 
vious conclusion about the disappearance of the small 
strained sized rings which react with the sodium ions 
during their adsorption at the amorphous silica surface. 
As observed for the ring size distributions, the orienta- 
tion of the rings (FIG.QJb)), in the three regions of the 
sodo-silicate film is similar to that obtained in the inte- 
rior of the pure silica films 2 ^ which means that even at 
the surface the rings have an isotropic orientation with 
respect to the surface. This is related to the disappear- 
ance of the small sized rings at the surface. Nevertheless, 
it is worth noting that in the external region the prob- 
ability of a Si atom to belong to a 5-fold ring is greater 
than the probability to belong to a 6-fold ring (FIG. 01(a)) 
as observed at the ASF surfaces 2 ® . In a sense, the 5-fold 
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FIG. 4: Ring probabilities (a) and mean orientations (b) as a 
function of their size for the • interior, ■ intermediate and O 
external regions. The same quantities calculated for the ASF 
surface ( — ) and ASF interior (■ ■ ■ ) regions—. 



rings are not affected by the sodium adsorption in con- 
trast to the small sized rings which disappear with the 
introduction of the sodium atoms. Also in the external 
region, the 2 and 3-fold rings (that are still present) are 
oriented perpendicularly to the surfaces (FIG. Efb)) as 
observed at the ASF surface d 28 ! 49 . 

We have also analyzed the dynamics of the Na + cations 
and compared the results obtained in the present 'NS10' 
system with those obtained in a NS4 glass. To this pur- 
pose, the mean square displacements (MSD) (r 2 (t)) = 
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\ri{t) — ?'i(0)| 2 y have been calculated for each species 

composing the SSF. FIG. |S] represents the MSD for the 
BO, NBO, Na and Si atoms within the time frame 0.7 fs - 
70 ps after the first 210 ps together with the MSD calcu- 
lated by Sunyer et al£& for NS4 {t ^ 1.5 ps). 

At ~ 2800 K, the MSD of each species exhibits three 
regimes. For the short times, we observe the so-called 
ballistic regime where (r 2 (i)) ~ t 2 . In this regime, the 
differences between the species are not really important. 
For the long times, we recognize the so-called diffusive 
regime in which (r 2 (i)) ~ t and in which the Si atoms 
are the ones that diffuse the less. Similarly to the pure 
ASF 28 , the NBOs diffuse more than the BOs. The 
origin of this feature lies in the fact that the NBOs form 
only one covalent bond with the Si atoms instead of two 
for the BOs. The Na atoms diffuse much more (one 
order of magnitude) than the other species. Between 
the two former regimes, the MSD exhibits the so-called 
/3-relaxation, clearly observable for the oxygen and 
silicon atoms. The phenomenon is not as clear for the 
sodium atoms, but between ~ 10 _1 and ~ 1 ps, the 
Na MSD does not behave like t (sub-diffusion) and it 
is reasonable to consider that the sodium atoms are 
submitted to the so-called cage effect characterizing the 
/5-relaxation. 




t[ps] 

FIG. 5: Mean square displacements for the BO ( — ), NBO 
(•••), Na (- -) and Si (• - ■) at ~ 2800 K. The thin lines 
represent the results obtained by Sunyer et al~— for a NS4 
system at 3000 K. 

These results are consistent with the observations 
made for different sodo-silicate systems^iSi* 2 ^ for tem- 
peratures close to 2800 K. More precisely the mean 
square displacements in the SSF are close to those cal- 



culated by Sunyer et alM for a quenched NS4 glass at 
3000 K. This is somewhat surprising since in a NS10 sys- 
tem, one expects smaller diffusion constants for all the 
species compared to the diffusion constants obtained in 
a NS4 system. The modus operandi of the present sys- 
tem (diffusion of the cations through the surface) could 
explain this observation. 



IV. CONCLUSION 

This work was motivated by a recent experimental 
study^i of the local environment of diffusing sodium 
atoms deposited at the surfaces of thin amorphous sil- 
ica films. We have reproduced numerically this experi- 
ment by classical molecular dynamics simulations after 
putting Na20 groups at the surfaces of amorphous silica 
thin films. 

We have quantitatively analyzed the temporal evo- 
lution of the sodium and total densities and we have 
checked that the sodium atoms are diffusing inside the 
amorphous silica network. After a given time, the den- 
sity profiles are no longer evolving, and we have calcu- 
lated the structure of the resulting sodo-silicate glass. 

Our attention has been focused on the local environ- 
ment of the sodium atoms. Once inside the thin film, 
they are preferentially bound to NBOs as predicted by 
the MRN model of Greaves^. The distances and bond 
angle distributions show that the sodium atoms have 
a local environment corresponding to the local environ- 
ment of the sodium atoms in sodo-silicate glasses made 
by quench, as observed by Mazzara et al, 2i . Moreover, 
the distances dNa-o = 2.2 A and rfsi~Na — 3.5 A are 
close to the experimental values. 

Concerning the amorphous silica network, we have 
observed via the corresponding distance and bond angle 
distributions that its short range order is not modified 
by the introduction of the Na + cations. We have also 
calculated the ring size distributions and the orientation 
of the rings which show that the introduction of the 
sodium atoms has an influence on the silica network 
but on larger scales compared to those corresponding to 
the local environment. The ring size distributions and 
the orientations are close to the results obtained in the 
bulk of thin amorphous silica films. This is due to the 
decrease of the proportion of small rings (particularly 
two and threefold) which interact with the Na20 groups 
since they are known to be highly reactive sites for the 
adsorption of species on amorphous silica surfaces. 
Finally concerning the dynamics of the different atoms 
we find results similar to those obtained in a NS4 glass 
at a slightly higher temperature. 
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